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THE EQUILIBRIUM OF A THIN COMPRESSIBLE MEMBRANE 
WITH APPLICATION TO THE PERIODONTAL MEMBRANE! 


By G. E. Hay? 


Abstract 


The equilibrium of a thin compressible elastic nnembrane has been considered. 
In the present work the theory developed in that paper (3) is applied to the perio- 
dontal membrane which lies between the tooth and the jaw-bone. 


For simplicity, a tooth is considered of which the root is a right circular cone, 
the thickness of the membrane being assumed constant. A load is applied to the 
tooth ina plane containing the axis of the root and a determination is made of the 
resulting stress in the membrane and displacement of the tooth. 


It is also found that the periodontal membrane holds the tooth in place in a 
way different from that commonly accepted by dentists. 


Introduction 


In a preceding paper (3) the equilibrium of a thin compressible membrane 
has been considered. In the present paper the theory there developed is 
applied to the periodontal membrane which lies between the tooth and the 


jaw-bone and is attached to both. The theory is applicable to multiple- 
rooted as well as single-rooted teeth but only the latter type will be con- 
sidered, in particular the upper central incisor which has the following average 
measurements (1, p. 19) :— 


Length over all 0.88 in. 
Crown length 0.39 in. 
Root length 0.49 in. 
Mesio-distal diameter at neck (side to side) 0.24 in. 
Labio-lingual diameter at neck (front to back) 0.27 in. 


The ayerage thickness of the periodontal membrane is from 0.0091 to 0.0099 
in. (6, p. 123). 


The equilibrium of the periodontal membrane considered as incompressible 
has been considered by J. L. Synge (8), and there is an interesting comparison 
between his results and those of this paper. 


In the author’s preceding paper (3, §3) there was introduced the equation 
(1.1) 1—20=ke"’, 
1 Original manuscript received September 28, 1938. 


Contribution from the Department of Applied Mathematics, University of Toronto, 
Toronto, Canada. 


2 Holder of a Bursary under the National Research Council of Canada. 
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where o was Poisson’s ratio, € was a small dimensionless parameter defined 
by 
(1.2) 


_ average thickness of the membrane 
~ average length of the membrane’ 


and m was an integer to be so chosen that k was comparable with unity. 
Owing to the method of approximation used in solving the equations of 
equilibrium, different results were obtained in the four cases »=0, 1, 2, and 
n>2. The first step is thus to determine which case applies to the perio- 
dontal membrane, for which the following average values may be taken 


e=2xX10-, e=4xX10-, e=8xX10-, 
(1.3) 

M=rigidity (2) =73 lb. in. 
Since the human body is composed largely of water it seems logical to take 
for the bulk modulus of the membrane that of water 


(1.4) K=2.96X105 Ib. 
Hence 

(1.5) 1 2.510, 


and so if m=1,2,3, the respective values of k are 


(1.6) k=(1—20)/e=0.6, 


k=(1—20)/e=30. 


Thus it appears that the case »=2 applies to the periodontal membrane. 
This case has been investigated (3, §8): it was found that the principal part 
of the stress in the membrane was a hydrostatic pressure satisfying a partial 
differential equation (3, Equation (8.22) ), subject to the boundary condition 
that the pressure at the free edge be atmospheric. This differential equation 
appears in a modified form in §2 and is solved in §3 for a restricted shape of 
tooth. 


2. The Tooth of Revolution. The Pressure in the Membrane 
Expressed by the Solution of Three Ordinary 
Differential Equations 


For a general tooth the partial differential equation satisfied by the pressure 
is rather complicated. The case will now be considered when the socket S 
and tooth T both have surfaces of revolution about a common axis. This 
implies that the thickness of the membrane has the same symmetry of revolu- 
tion. The edge of the membrane is assumed to be a circle normal to the 
axis of revolution, although in practice this is not realized very closely. 
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Fic. 1. Model of a tooth and periodontal membrane. 


A set of rectangular cartesian axes are taken with the origin at the apex, 
the axis of x being coincident with the axis of revolution. A set of surface 
co-ordinates ¢ and @ are also taken over S where ¢ is the meridian distance 
from the apex and @ is the azimuth measured from the plane z=0. Then 
the partial differential equation for the pressure p in the membrane [3, Equa- 
tion (8.22)] takes the form 


(2.1) (xi 32) +2 ~ 12(1—20)hp =4E?; 


the boundary conditions are the following, 
j 


for ¢=0, p must be single valued, 
(2.2) 
for p=P; 


and R is the radius of a general section normal to the axis of x, h is the thick- 
| ness of the membrane, E is Young’s modulus, ¢ is Poisson’s ratio, 6° is the 
component normal to S of the displacement of a point on 7, counted positive 
when in the sense from S to T, / is the meridian length of the membrane and 
P is atmospheric pressure. 


In (2.1) B° is a function of ¢ and @ and is known when the infinitesimal 
rigid body displacement of T is known. To find it, the displacement of T 
is resolved into infinitesimal translations u, v and w in the directions of the 
axes and infinitesimal rotations @;, Wz, and ws; about them. Then the dis- 
placement of a point (x, y, z) on T has components : 


U—YWs+2We , W—XWet 
and hence 
(2.3) dt cos { + ws (« dt +R 
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Thus (2.1) can be written in the form 


Ge AE cos — ws (x ae +R 


—4E sin w — ws (« +R 


It is to be noted that w; does not appear in (2.4), since it contributes nothing 
to B°. 

Except for cos @ and sin @ the right-hand side of (2.4) is a known function 
of ¢ and the solution of (2.4) must then have the form 


(2.5) b=f(S)+ cos sin 


Upon substitution of this in (2.4) a differential equation is obtained for each 
of f, g and g’, whence follows the result :— 


I. In the case of a tooth of revolution the pressure in the periodontal mem- 
brane is given by (2.5) where f, g and g’ satisfy the equations 


(2.7) - { Fe +12(1-20)h | 

(2.8) RE (Re +12(1—20)h } 


with the boundary conditions 


(2.9) f(0) must be finite, f(/) =P, 
(2.10) g(0) =g(/) =0, 
(2.11) g’(0) =g’(1) =0. 


In these formulae R, ¢, $, x and / are as shown in Fig. 1, h is the thickness of 
the membrane, E is Young’s modulus, @ is Poisson’s ratio and the infinitesimal 
rigid body displacement is specified by the translation (u, v, w) and the rotation 
(Wi, W2, Ws) relative to the axes of x, y and z in Fig. 1. 


| 
| 
a 
t 
| 
} 
] 
| 
| 
- | 
| 
Fe | 
: 
| 
| 
| 
f 
i 
an 4 


HAY: EQUILIBRIUM OF THE PERIODONTAL MEMBRANE 127 


3. The Conical Model under Plane Displacement 


A tooth of which the root is a right circular cone will now be considered. 
This is the simplest geometrical figure which can be chosen to represent the 
average shape of a single-rooted tooth. The rest of the paper is confined to 
membranes of uniform thickness 4, although the thickness varies consider- 
ably in the human tooth. 

Only plane displacements of the tooth will be considered, 7.e., displace- 
ments in which the axis of the tooth remains in a plane, and the axis of y 
will be chosen to lie in this plane. Hence in the notation of §2, 


(3.1) W=0,;=02=0. 
v} 
\ 
— > 4 x 
Fic. 2. Conical model of a tooth. 
Fig. 2 shows the conical model for which 
(3.2) sin a, cos a, 
and substitution from (3.1) and (3.2) in (2.6), (2.7) and (2.8) gives 
(3.3) + —12(1—20)hf =4Eu sin a, 
rg 
= —4E(v cos a+ wf), 
(3.5) de? dt {120 -20)h+ g’=0, 


A variable £1, a constant J, and the functions f:, gi and g,’ are defined as 
follows :— 


(3.6) 1,=vl, 
(3.7) =f(0), =g($), gr ($1) =2'(8), 
where 


(3.8) 
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Then the differential equations (3.3), (3.4) and (3.5) take the forms 


(3.9) 


(3.10) 


(3.11) 


The boundary conditions obtained from (2.9), (2.10) and (2.11) are 


(3.12) f:(0) must be finite, f:(/:) =P, 

(3.13) gi(0) =0, 

(3.14) g1'(0) =0, 

where G, C, and D are the constants 

(3.15) G=4Eu sin a/(h*v?), C=—4Ev cos a/(h*v’), = —4Ew;/(h*v*) , 


and all other quantities are defined in I. 


Solutions for (3.9), (3.10) and (3.11) will be obtained in turn. Since (3.9) 
is Bessel’s modified equation of order zero, the solution satisfying the boundary 
condition (3.12) is 


(3.16) 
where 


(3.17) 


In the course of solving (3.10) it is found that the form of the solution is 
different in the three cases (i) cosec a ~ integer, (ii) cosec a = odd integer, 
(iii) cosec @ 
arithmetical values in §1 give cosec a = 4.0. Thus a solution will be required 
in Case (iii) only, in which case 
(3.18) 


The left-hand side of (3.10) is then Bessel’s modified equation of order 2n 
and the complete solution of (3.10) is 


(3.19) 


where A and B are arbitrary constants, D and C are the known constants 
defined in (3.15), Ion(€1) and Kena({1) are Bessel functions and W and x are 
any solutions of 


(3.20) 


(3.21) 
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+ $1 —£2(G+f1) =0, 


+o, — (cosec?a+ =0. 


Gs)" 


= even integer. Since @ is the semi-angle of the root cone, the — 


cosec a=2n, 


=DW(S1) + CX ($1) +A + BK , 
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When n is an integer (3.20) has the solution in series* 


(3.22) Wh) = 
{(2r+1)?—4n?} 


and (3.21) has the solution in series** 


(3.23) 
II {(2r)*—4n} 


An {(2r)?—4n?} { (2r)?—4n?} 


© 1 m 4r 

r=1 


where the dash on the II and 2 signs means that the terms r=m are to be 
omitted. Using the boundary conditions (3.13), a determination can be made 
of A and B, whence 


(3.24) = + Cx ($1) — + Cx (Ar) } Len ($1) (21) « 
It is noted that 


Since (3.11) is Bessel’s modified equation of order 2n, it has the complete 
solution 
(3.26) 81'($1) =A’Tan($1) +B’ 


where A’ and B’ are arbitrary constants. Using the boundary conditions 
(3.14), it is found that A’=B’=0, whence 


(3.27) =0. 
These results are summed up as follows:— 
II. For a conical model of a tooth subjected to the plane displacement 


specified by u, v, and ws as shown in Fig. 2, the pressure in the periodontal 
membrane is given by 


(3.28) =fi(F1) + cos , 
© This is derived in Appendix A. 
** This is derived in Appendix B. 
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where 


= , 
= DW ($1) + Cx ($1) — {DW(h) + Cx } 
G=4Eu sin a/(h*v’), C= —4Ev cos a/(h*v*®), D=—4Ews/(hv*) , 


v=V12(1—20)/h, cosec a@=2n, (m=integer), 


and wW(f1), x(f1) and J2,(f1) are defined in (3.22), (3.23) and (3.25) 
respectively. 


In these formulae ¢ is the azimuth measured from the plane z=0 of the dis- 
placement (Fig. 2), P is atmospheric pressure, E is Young’s modulus, ¢ is Pois- 
son’s ratio, h is the thickness of the membrane and is a constant, a@ is the semi- 
angle of the root cone and ¢; is a co-ordinate running along a generator of the 
root cone, having the value zero at the apex and at the edge of the membrane 
the value /,;=vl, where / is the slant length of the root cone. 


The external force system, exclusive of atmospheric pressure, acting on 
the tooth and producing such a pressure » in the membrane is statically 
equivalent to forces X and Y acting along the axes of x and y respectively and 
a couple N acting about the axis of z: thus 


1 
me | (p—P)f sinter df dd , 


o=0 

(3.29) Y= (p—P)f sin acos acos df dd, 
t=-0 J¢=0 


N= (p—P)f? sin a cos df dd. 


o=0 


Since [:=v¢ substitution for p from II and easy integration gives the fol- 
lowing result :-— 


III. The external force system necessary to produce the pressure distri- 
bution in II is statically equivalent to the force system (X, Y, N), (Fig. 2), and 


(3.30) X= (P+G) sin® 
(3.31) ¥= sin acos a{Dya(h)-+Cxa(h) —[D Wh) + Cx(h) 


(3.32) N= sin a{ + Cx2(h) —[DW(h) + Cx(h) en2(h)/Ten(h) 
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where 


ly 
(3.33) Teon,t(h) = dh, 
0 


=2tH 
m=0 (2m+2n+t+1)m!(m+2n)! ’ 


(O!=1, m=0,1,..., #=1, 2), 


ly 
(3.34) -[ dh, 
0 


],2™+ t+2 


=> 
™=! (2m+t+2) {(2r+1)?—4n?} 


’ (¢=1, 2), 


(3.35) = i dh, 
0 


= t+ 
™=1 (2m+t+1) {(2r)?—4n?} 
2n+t+1 2m+t+1 
4n {(2r)?—4n?} (2m+t+1) Th {(2r)?—4n*} 


and the notation is as in II, the dash on the 2 and II signs again indicating 
that the terms r= are to be omitted. 


The following values consistent with one another and with the physical 
data are assumed as the basis for calculation: 


b=over-all length of tooth =0.88 in. 
1=slant length of root cone =0.50 in. 
(3.36) a=semi-angle of root cone = 
h=thickness of membrane =0.0098 in. 
v=+V/12(1—20)/h =5.6 in. 


In §5, on account of the cancellation of positive and negative terms, it is 
found necessary to work initially with eight-figure accuracy in order to obtain 
three significant figures in the result. 


q 

. 
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4. The Conical Model Subjected to a Force X Acting on the 
Crown and Parallel to the Axis of the Root 


The crown is assumed to lie on the axis of the root. The tooth then under- 
goes an axial displacement only, and in the notation of §3, 


(4.1) Y=N=0, v= C=D=0. 
Since G=4Eu sin a/(h3v*) the axial displacement is, from (3.30), 


Xv? 
and in terms of u the pressure distribution is, from II, 
_  4Eusina 4Eu sin a\ Io(f1) 


It is to be noted that the membrane is in the unstressed state when the 
external load and the atmospheric pressure both vanish. When X=0 the 
whole system (tooth, membrane and jaw) is subject to atmospheric pressure 
P only, and the axial displacement, which will be denoted by %, is 


(4.4) —Ph'y?/(4E sin a). 
Then from (4.3) and (4.4), 
(4.5) (£1) = sin =P, 


i.e., the pressure throughout the membrane is atmospheric. 

The points of maximum pressure in the membrane will now be determined 
for a general axial load on the tooth. The pressure at the edge of the mem- 
brane is always atmospheric, and, since the gradient of Jo({1) is positive, 
from (4.3) it can be seen that 


4 >0 or <0 according as P+4Eu sin a/(h*v®) >0 or <0. 


Substitution for u from (4.2) gives 
>0 or <0 according as X } >0 or <0. 


Using the values in (3.36), it is found that J;=2.80, =4.16 and Io = 
9 24, whence /To(l1) =3, 40. Thus 


3 >0 or <0 according as X>0 or <0. 


These results are summed up as follows:— 


IV. When the external load acts from apex to crown the minimum pressure 
occurs at the apex and the maximum at the edge of the membrane. When 
the external load acts from crown to apex the minimum pressure occurs at 
the edge of the membrane and the maximum at the apex. In each case the 


& 

* 
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pressure at the edge of the membrane is atmospheric and the pressure at the 
apex is 
Xv? 


sin? a 


(4.6)  p(0)=P-— 


{1—1/Io(h)} , 


where P is atmospheric pressure, v=5.6 in.—!, 1;=2.8, Io(li) =4.16, Io, 
=9.24 and sin a=}. 


A diagram will now be constructed illustrating the pressure variation along 
a generator of the root cone when X is such that, (7) p(0) =2P, (4) p(0) =0. 


When (0) =2P, the external load will be denoted by X’, the corresponding 
displacement by u’ and the corresponding pressure by #’(f{1). From (4.6) 
it is found that 


(4.7) X'=- sin?a {l,?— 2Ioa(h)/Io(h)} 


and hence by (4.2), 


(4.8) “=~ GEsina Ioh)—1 


Finally from (4.3), 
(4.9) (£1) =P {2Io(11) —1 — } { —1} +. 


When #(0) =0, the external load will be denoted by X”’, the corresponding 
displacement by u” and the corresponding pressure by p’’({1). It is found 
that 


(4.10) X"=—X’, b''($1) =2P—p'($1). 

To compute u’ and u’’ requires Young’s modulus for the membrane. Since 
o =}, approximately, then E=3y: 
since w is known (2) 
(4.11) E=219 Ib. in. 
When P= 15 Ib. in.~ the values in (3.36) give 

X'=—X" =—0.42 lb., 

(4.12) = —u"’ = —2.65X10-* in., 


a= —1.99X10-* in. 


The graphs of p’({1) and p’’({1) are shown in Fig. 3. 
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When the membrane is incompressible the corresponding values derived by 
Synge (8) are 


X’= —-X" = —0. 38 Ib., 
(4.13) 
= —u'’=—8.8X10-7 in., 


and the graphs corresponding to p’({1) and p’’({1) are shown in Fig. 3. 


2P Compressible | Incompressible 
membrane membrane 
B; 
! x -0.42 lb. -0.38 lb. 
-8.8x1077in. 
Pressure 
A B 
dist'n. ' 
uJ 
0.00 lb. 
Ag -1.99 x107%1n, 
uJ 
a Pressure A 
xX 0.42 1b. 0.38 1b. 
u 6.6x107 in, 8.8x 107 "in. 
Pressure 
Bz . dist'n. 


P =atmospheric pressure (=15 Ib. in.~*), 
X =external load, 
u=displacement due to X and atmospheric pressure. 


In the case of the compressible membrane the total displacement u is made up 
of a part due to atmospheric pressure and a part due to the load X. When 
X =—0.42 lb., the part due to atmospheric pressure is —1.99X10~ in. and 
the part due to X is —2.65 107 in. 


Fic. 3. Conical model under axial load. 


5. The Conical Model Subjected to a Force Y Acting on the 
Crown and Perpendicular to the Axis of the Root 


The calculations in this section must be carried out initially to eight figures, 
as explained at the end of §3. The crown is again assumed to lie on the axis 


xX 
a 
- § 
x 
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of the tooth. Then, in the notation of §3, X =0 and N=DY, where 3 is the 
over-all length of the tooth. Hence from (3.30), (3.31) and (3.32), 


(5.2) Cai+DB:= sin cos 
(5.3) Caz+ DB.=bYv*/(m sin a), 
where 


= Wall) — Tan, /Ten (li) , 
5.4 


= X2(L1) —X (11) — , 
and the remaining notation is defined in II and (3.33), (3.34) and (3.35). 


Numerical values are now introduced from (3.36). Then n=2 and the 
following values can be computed :— 


( 


=0.234 079 08, (li) =— 6.601 251 5, 
x(l1) =—1.618 780 30, 
=0.279 338 86, villi) = 8.601 708 7, 


(5.5) =—2.309 403 8, 
=0.677 170 52, = —20.547 149, 
=—5.407 079 8, 
whence 
—0.377 628 3, Bi=—0.724 088 6, 
(5.6) 
O2=—0.724 088 5,  B2=—1.450 300. 


The solution of (5.2) and (5.3) is thus 
(5.7) C/Y=3 537.441, D/Y=~—1 901.81, 
and using (5.1), the pressure distribution is, from II, 


(5.8) “p(t:)—P=Ycos ¢[—1901.81 (61) +3537 .44x (61) —2917014(¢,)], 


where ¢, runs from zero at the apex to 2.8 at the free edge, @ is the azimuth 
as shown in Fig. 2 and ¥({1),x(1) and J4({1) are defined in (3.22), (3.23) and 
(3.25) respectively. The graph of p({1) is shown in Fig. 4. The maximum 
pressure occurs at M and N, and it is found that the external load necessary 
to produce a maximum pressure of two atmospheres (30 Ib. in.~*) is 


(5.9) Y=0. 188 lb. 


From Synge (8) the corresponding load for the incompressible membrane is 


Y=0.190 lb. 


(5.1) | P+G=0, 
. 
‘ 


2P, 


r C,fCe2 \ 


PRESSURE 


P =atmospheric pressure (=15 Ib. in.~*), 
Y =external load, 


v=displacement produced by Y. 
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Compressible Incompressible 
membrane membrane 
Apex to centre of rotation 0.29 in. 0.31 in. 
Zz 0.188 Ib. 0.190 Ib. 
v 3.62 X10- in. 2.8107 in. 
Pressure distribution on mn CQ C2 


Fic. 4. Conical model under transverse load. 


The definitions of C and D are given in II, whence from (5.7), 


(5.10) v= —35.4X10? 


Vhiy? 
4E cosa’ 


w3=19.0X10? 


4E 


where the units of measurement are the pound and inch. 
Now £=219 Ib. in.~*, and when the biting edge is acted upon by the critical 


transverse load (0.188 Ib.) it undergoes a displacement 6 and 


(5.11) 


’ 


=3.62X10-5 in. 


The corresponding displacement in the case of the incompressible membrane 


is 6=2.80X10- in. 


| 
Were 
\ 
\ 
n 
Vv 
; 
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The centre of rotation of the tooth is at a distance x’ from the apex where 
v+x’w;=0, whence x’ = —v/ws; and from (3.36) and (5.10) it is found that 
x’=0.29 in. Since the axial length of the root is 0.48 in. the following result 
is obtained. | 


V. When the conical model is subjected to a transverse load applied at 
the biting edge, the centre of rotation divides the root approximately in the 
ratio of two to three and is farthest from the apex. 


This result agrees roughly with that of R. Kronfeld (5, p. 1270). 


6. Discussion 


There are blood capillaries in the periodontal membrane. The blood 
is squeezed out of those parts of the membrane where the blood pressure is 
exceeded by the pressure due to the external load on the tooth. There is 
said to be strangulation of the membrane in such regions. From Fig. 3 it 
is seen that, if axial load is steadily increased from zero, strangulation will 
first appear at the apex of the root. From Fig. 4 it is seen that, if transverse 
load is steadily increased from zero, strangulation will first appear at M and N. 
Computation shows that strangulation will be produced in small regions 
surrounding these points by a transverse load of 3 gm. Schwarz (7) con- 
ducted experiments on the teeth of dogs and found that the symptoms of 
strangulation appeared for a transverse load of 5 gm. This agreement is 
quite good, as the teeth of dogs differ in shape from those of humans. 

The periodontal membrane contains a large number of small fibres running 
from the tooth to the jaw-bone. It is commonly believed among dentists 
that the resistance against displacement of a loaded tooth is supplied by 
tension in these fibres (4, pp. 336-337). In this paper it is found that the re- 
sistance is attributed to the pressure in the membrane, and the agreement 
in the previous paragraph between calculated and empirical results substan- 
tiates this claim. 
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Appendix A : 
To obtain a solution in series of the equation 
2 | 


df? 


where n is an integer. 
This equation has a solution of the form 


(A.2) y= z 


- 
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if b, satisfy the equations 

b3(3?—4n?) =1, 
(A.3) b4(4?—4n?) =0, 


If 

(A.4) (s=4,6,...); 
then a solution of (A.1) is 

(A.5) y= 


m=! { (2r-+1)?—4n?} 


Appendix B 


To obtain a solution in series of the equation 


where n is an integer. 
If an attempt is made as in Appendix A to obtain a solution of the form 


(B.2) 2 
it is found that | 


(B.3) = ©, (¢=0, 1, 2, 


and so some modification is necessary. 
The equation 


2: 
(B.4) +6, — e<<i, 


has a solution of the form 


if 
c2{(2+€)?—4n?} =1, 


(B.6) c3{(3+6€)?—4n?} =0, 

—c,2.=0, Gee .« 
Thus 
(B.7) c,=0, (eed, 


4 
© 
4 
te 
~ i 4 
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and so 


{(2r+)?—4n?} 


Since (2n+ €)?—4n?=e(4n+e), then 


IL 

where 

(B.10) 


(4n+e) { (2r+€)?—4n?} 4n?} 
Now 


whence 
(B.12) 


r=1 


Similarly by changing the sign of € a function x. can be obtained. Hencea 
solution of 


is 


(B.14) 
{(2r+)?—4n?} {(2r—€) —4n?} 
+$(0) log. tel... 


By letting e— 0 simultaneously in (B.13) and (B.14) it is found that a solution 
of (B.1) is 


X= lim. x, 
e—0 


n—1 62" 
- 5 +$(0) loge 
m=1 { (2r)?—4n?} (3 
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or 


(B.15) 
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log. 


4n i {(2r)?—4n?} 


meat! {(2r)?—4n?} 
r=n+1 


cr 


Laat 2 ] 


m=n 4n r=1 (2r)?—4n? 


where the dash on the Z andII signs means that the terms r = 7 are to be omitted. 
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THE RESONANT RADIAL FREQUENCIES OF A CYLINDER 
OF ANY WALL THICKNESS! 


By G. S. FIE.Lp? 


Abstract 


Equations have been derived for the resonance radial frequencies of a cylinder 
of any wall thickness. For the degenerate case of the internal radius equal to 
zero, the equations are shown to reduce to those previously given for a solid rod. 
For a cylinder with a very thin wall, the frequency equation becomes very simple, 
but contains an additional factor compared with the equation usually given fora 
thin-walled cylinder. It is suggested that the two cases correspond to slightly 
different types of vibration, and that the one considered in this paper is physically 
more likely. 


Radial Particle Velocity 


ies a previous paper (1) it was shown that for an element of a solid cylindrical 
rod, the expression for the radial particle velocity, is, 


= , 


where, ¢ = radial particle velocity at a distance r from the centre of the rod; 
B = constant of integration; 
w = 2m X frequency; 
t = time; 
z = distance along the length of the rod; 
¢: = phase velocity of longitudinal wave in the rod; 

Ji(kr) = Bessel function of the first kind, and is a solution of the following 
equation: 
2 
— 1)fe = 0, 


in which k is a function of c,, w, Young’s modulus (g), the density (p), and 
Poisson’s ratio (7), but is independent of r. 

The complete solution of Bessel’s equation, above, includes Bessel’s function 
of the second kind, 7.e., Yi(kr)*, but for the solid rod this solution is not 
admissible, as Yi(kr) becomes infinite at r = 0. 

For a hollow cylinder, however, Yi(kr) must be included, and the expression 
for the radial particle velocity becomes, 


= [AJi(kr) + (1) 


Axial Force on an Element 


The equation for the axial force, 2, along a radius vector, on an element 
in the wall of the cylinder, may be put in the following form (1, p. 257), 


1 Manuscript received March 28, 1939. 
Contribution from the Division of Physics, National Research Laboratories, Ottawa, 
Canada. Issued as N.R.C. No. 821. 
2 Physicist, National Research Laboratories, Ottawa. 


* Y,(x) was adopted as the canonical function of the second kind by Weber (3, p. 64) and is 
the most useful in the present case. 


c 
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where, 
& = longitudinal particle velocity at a distance z along the axis; 
n = tangential particle velocity corresponding to the cylindrical 
co-ordinate 0. 
If only radial vibrations are considered, —- > 0, and, 


Boundary Conditions 


Let the external radius of the hollow cylinder be a; and the internal radius 
be a2. Then, 


Atr = 0; 
Atr = a2,f2 = 0. 
Substitute in Equation (3) for f andr. Then, 


Alkr(1 — o)Jo(kr) — (1 — 20)Ji(kr)] 


= 0. (4) 
+Blkr(1 — o)Vo(kr) — (1 — 20) Vilkr)] 
Let kay; = «x, 
Then kaz = (a2/a;)x. 
Divide each equation by B(1 — ¢@). 
Write,m and = A/B. 
Finally, 
FlxJo(x) — mJi(x)] + xYo(x) — mYi(x) = 0, (5) 
and 


— mJi(a2x/a)] 
+ (d2/a1)xVo(aex/a1) — mYi(aex/ai) = 0. (6) 
From these equations it is possible to determine F and x. 


Frequency Equation 
The expression determining the resonant radial frequencies is (1, p. 259), 
f = ka/2r, (7) 
frequency of resonant radial vibration; 


where, f 


modulus of elasticity 
density 


gi — 
(1+ — 


3 
may 
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Since ka; = x, 
xa 
f (8) 


Solution of Equations in x and F 
In Equations (5) and (6), F may readily be eliminated; this leaves an 
equation containing only x and functions thereof. It is possible to solve 
this equation for any ratio of a2/a; by substituting for different values of x 
the values of the Bessel functions as given in the published tables. 
A simpler method of solution is a graphical one. In Equation (5) denote 
F by F;, and in Equation (6) write F as Fe. 


_ _ — mYi(x) _ _wN 
xJo(x) — mJi(x) (9) 
— mJi(a2x/ai) G 


From these equations, values of x may be determined from the intersections 
of F; and Fs. By a consideration of the change-over points, where either 
the numerator or denominator of each equation equals zero, that is, where 
F; or Fe changes sign, the ranges of x where solutions are possible may be 
very greatly reduced. As an illustration of the method, definite values of 
m and of the ratio a2/a; may be taken. 


Consider the values of x for which solutions are possible. 
Assume 0 = 0.226, from which m = 0.709, and assume az/a; = 0.8. 
Table I may then be set up. 


TABLE I 
Possible 
solution 
x N D H G F; Fs (Fs a Fe) 
marked + 
0 - 2.04" + + + + = = + 
2.04 2.55 + 4 4: 
2.55 + 0 + —© 
3.76 — 4.70 + + 
. 4.70 0 0 
etc. etc. etc. etc. 


From this table it will be observed that there are only a few ranges of x 
where solutions are possible. It would be expected that in the first range 
would be found the solution giving the fundamental resonance frequency; 


. 
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in the second, the solution for the first overtone; in the third, the second 
overtone; and so on. 

By plotting a few points for F; and F, in each range, the form of each curve 
may be readily obtained. From the intersection of the curves, the desired 
solution of x is determined and hence the resonant frequency. In Fig. 1 
the curves for a particular case are given. 


1-2 2: 
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INSIDE RADIUS (a,) 
Fic. 1. Curves of Fsand Fe , illustrating graphical solution for x. =0.709, a2/a; = 0.8. 


Fic. 2. Curve showing changes in the resonant radial frequency of a cylinder with changes 
in the internal radius. Cylinder made of —_ external radius, 1.28 cm. Constants of glass 
as follows: g = 6.4 X 10", p = 2.6,0 = 0.226. 


A number of resonant radial frequencies were calculated for values of x 
determined as above, to show the variation in frequency with variations in 
the ratio a2/a,. In Fig. 2 these are plotted against az. 


Limiting Case of the Solid Rod (a, — 0) 
In Fe write (a2/a;)x = y. 
When az 0, y > 0, Jo(y) 1, ~ 0 and 


: 
9 
| 
. 
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Hence Fs — ©, unless H — 0, in which case Fs, becomes indeterminate. 
Suppose H > O0asy— 0. Then, 


When y > 0, Yo(y) ~ —. Hence, for y — 0, yYo(y) can be put in the 
form o/o, Remembering that for this indeterminate form the limit of 
the quotient is equal to the limit of the quotient of the derivative of the 
numerator and the derivative of the denominator, the following may be 
obtained. 


Then, 


But, since the limit of the product of a finite number of quantities is equal 
to the product of their limits, 


Hence, 


Lim _ Lim 
y — 0 (™) y 


which is obviously untrue. 
Therefore, as ag — 0, H does not — 0, and Fy > ~. 
Since F; = Fs, asaz—0, Fs > 
Hence, for a. 0, 


xJo(x) — mJi(x) 0, 
xYo(x) — mY,(x) @. 


Since xYo(x) — mY,(x) can approach infinity only for x — 0, the only 
useful solution for x is contained in the expression, 
xJo(x) mJi(x) . 


It has previously been shown (1, p.259) that the frequencies of radial 
resonance of a solid rod are given by the equation, 


xJo(x) = mJi(x), 


so that for the limiting case (a2 — 0) the equations developed in this paper 
are seen to degenerate to the correct form. 


Lim 
y—0 
° 
Lim 
y—0 
Lim 
y—0 
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Limiting Case of the Cylinder with Infinitely Thin Walls (a, — a;) 
Write (a2/a:)x = x — h. 
Then, h = (1 — a2/a:)x, 

and, as > 0. 


From Taylor’s theorem, 
h? 
— h) = fx) — @) + 
By means of this equation, Fs becomes, 


(x — A)[Vo(x) — hYo'(x)] — ml Vi(x) — hVi'(x)] 
(x — h)[Jo(x) — hJo'(x)] — miSi(x) — hJi'(x)] 


in which expression h? and higher powers of h have been neglected, since only 
those values of # that approach zero are considered. 


To obtain the desired solution for x, Fs is to be equated to F;, as follows. 
To simplify writing the equations, Jo is written for Jo(x), Yo is written for 
Y,(x), and so on. 


(x h)[ Yo hY'| mY; xVo mY, (11) 
(x = h)[Jo = hJ,'| xJo — mJ," 


In the preceding equation, if the various derivatives are written in and the 
terms in Jo, J:, Yo, and Yj, are collected, the following is obtained. 


mh 
+ mh) Yo + (hx hk? — m _ mY, 
This equation reduces to, 
| + mx —h+ mi) 
mh 
= Jo¥s{ x(x — — m — + — mi} (13) 
Hence, 
JoNi, (14) 
or, 
— — m — 2) + — mh) = 0. (15) 
The only solution of Equation (14) is x = ©, because, 


¥n(x) — Ja(x) Yn-1(x) = 2 


wx” 


Hence only Equation (15) need be considered. 


< 
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Equation (15) reduces to, 
x? = m(2 — m), (16) 
if h? is neglected. 


The radial frequency is, 


— — a) 1 — 20 _1-2 
since, 


Equation (17) reduces to, 


(18) 


This is not quite the equation ordinarily given (2, p. 138; for the case s = 0) 
for the frequency of free radial vibration of a thin-walled cylinder, as it con- 
tains the additional factor, Vea . As the value of Poisson’s ratio usually 
varies between } and 3, the extra factor leads to an increased frequency of 
the order of about 3 to 6%. Practically this is of no very great significance, 


but the theoretical reason for the difference is interesting. 


In the derivation of the above equations, a wave going from the inside of the 
cylinder wall to the outside of the wall, thence being reflected back towards the 
inside, and then being further reflected, was considered. In the derivation 
of the equation for the thin-walled cylinder, a simple extensional vibration 
of the tube wall is assumed; the phase of any element in the wall along the 
same radius vector is the same. 

From what has just been said it will be realized that the type of vibration 
is essentially different in the two cases. Physically, the radial vibration of a 
cylinder must always be excited by an alternating pressure on the inside or 
on the outside of the wall, so that it is the type of vibration discussed in this 
paper that must occur. The frequency given by the simple theory is thus only 
an approximation to the correct value. 
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A METHOD FOR THE DETERMINATION OF ORGANIC 
PHOSPHORUS IN SOILS AND SOIL EXTRACTS' 


By C. L. WRENSHALL? AND W. J. DyeEr® 


Abstract 


Dissolved organic phosphorus may be determined, even in highly coloured 
alkaline soil extracts, by the application of photoelectric colorimetry to the 
ceruleomolybdate colour reaction of phosphate ions. Total organic phosphorus 
may be extracted from soil with 4 N hydrochloric acid followed by treatment of 
the residue with 3 N ammonium hydroxide. Inorganic phosphorus is determined 
by direct colorimetry in the acid extract and in the acidified alkaline extract. 
Total phosphorus is determined similarly in each extract after ignition. The 
difference between total and inorganic phosphorus is equivalent to the organically 
combined phosphorus content of the soil sample. 

The soil organic phosphorus was found to be stable under these treatments. 
A considerable proportion of the organic phosphorus is extracted by acid, and the 
total soil phosphorus is almost completely extracted by the two treatments with 
4 N hydrochloric acid and 3 N ammonium hydroxide. 


Introduction 


One of the recent trends in soil chemistry is to resume the study of the soil 
phosphorus that is combined in organic compounds. This subject is being 
approached with two objectives; first, to obtain a better knowledge of the 
chemical combinations in which the soil phosphorus exists, and second, to 
explain the role of organic phosphorus in the complex cycle of phosphate in 
the soil. A review of the literature shows that no completely satisfactory 
method for the determination of organic phosphorus in soils and soil extracts 
has been proposed. 


Parker and Fudge (6) applied the ceruleomolybdate colour reaction of 
phosphate ions to the determination of organic phosphorus in water extracts 
of soils.. Inorganic phosphate was determined directly in an aliquot, and 
total phosphorus in another aliquot after ignition. The difference was taken 
to be organic phosphorus. This method was very limited in its application, 
since for visual colorimetry it was necessary to have colourless solutions. 

Recently Dean (2) employed the same principle in determining organic 
phosphorus in alkaline soil extracts which he decolorized with kieselguhr 
or by treatment with bromine. Experiments in this laboratory indicate 
that complete decolorization is not attained by the action of bromine, and 


1 Manuscript received February 16, 1939. 
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Quebec, Canada, with financial assistance from the National Research Council of Conale. 
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that the treatment with alkali and bromine is sufficiently drastic to release 
some phosphorus from organic combination. 


Various methods have been proposed for determining total soil organic 
phosphorus by using the ceruleomolybdate reaction in conjunction with 
treatments designed to release phosphate from its organic linkage. Odynsky (5) 
determined the increase in acid soluble phosphate produced by igniting the 
sample. Recently Dickman and De Turk (3) published the details of a 
method based on releasing the organic phosphorus by treatment with hydrogen 
peroxide. In both these methods the whole soil is subjected to the destructive 
treatment and there is a danger that released phosphate may become fixed 
by the soil residue. Davis (1) has demonstrated that some soils may absorb 
phosphate from a medium with a pH value as low as 0.25. The usefulness of 
these methods is limited because they may not be applied to coloured solutions, 
such as alkaline soil extracts, where direct visual colorimetry is impracticable. 


In a previous paper (4) the authors showed that it is possible to differentiate 
between organic and inorganic phosphorus in various soil extracts, including 
those containing coloured organic matter, by the use of a photoelectric colori- 
meter. It was thought that if all the organic phosphorus of soil could be 
brought into solution by suitable means, a determination of the total organic 
phosphorus could be effected. Such an approach proved feasible, and resulted 
in the method described herein. The method has been used with practical 
success since September, 1937. 


A. Extraction Procedure 


A 1 gm. sample of finely ground soil (100 mesh) is placed in a 50 ml. Erlen- 
meyer flask and treated with 25 ml. of 4 N hydrochloric acid (redistilled). 
The contents are mixed by gentle rotation several times during the first 
hour, and then allowed to stand overnight. The mixture is filtered through 
hardened filter paper, and the residue washed with hot 0.01 N hydrochloric 
acid, it being completely transferred to the filter in the process. Washing 
with pure water sometimes causes peptization of organic matter. The filtrate 
is collected in a 100 ml. volumetric flask, cooled, and diluted to volume. 

The residue from the acid extraction is then washed into a 100 ml. volu- 
metric flask containing 20 mi. of concentrated ammonium hydroxide, and 
diluted to volume. The flask is stoppered and the contents are mixed several 
times during an hour, and allowed to stand until the following day. A quan- 
tity of this extract is centrifuged at moderate speed (3000 r.p.m.) for a few 
minutes to remove suspended material before aliquots are withdrawn for 
analysis. 


B. Preparation for Total Phosphorus Determination 

A 10 ml. aliquot of the extract is placed in a silica evaporating dish con- 
taining 10 ml. of 1 N magnesium nitrate solution, and evaporated to dryness. 
Organic matter is completely destroyed by ignition at 550°C. After the dish 
has cooled sufficiently, 10 ml. of 1 N hydrochloric acid is added, and the 
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covered dish is placed on the steam bath for two hours to bring about com- 
plete hydrolysis to orthophosphate. The contents of the dish are then washed 
into a 100 ml. volumetric flask and one drop of paranitrophenol indicator 
(0.5% solution in water) is added. The pH of the solution is adjusted to 
approximately 3.0 by adding dropwise 1 N ammonium hydroxide until the 
appearance of a yellow colour, then discharging the yellow colour with a 
drop of 1 N hydrochloric acid. The solution is diluted to volume. A suit- 
able aliquot, usually 25 ml., is transferred to a 50 ml. volumetric flask, and 
2 ml. of the Truog and Meyer (10) acid molybdate reagent is added; this 
final solution is diluted to volume, and is then ready for the colorimetric 
determination. 


C. Preparation for Inorganic Phosphorus Determination 


(1) On the acid extract. A small aliquot, usually 5 ml., of the acid extract, 
is pipetted into a 50 ml. volumetric flask, diluted to 30 to 40 ml., the pH 
adjusted to approximately 3.0, and finally diluted to volume after the addition 
of 2 ml. of acid molybdate reagent, as described above. 

(2) On the alkaline extract. A 10 ml. aliquot of the centrifuged extract is 
acidified by the dropwise addition of 4 N hydrochloric acid with continuous 
stirring. Acidification is complete when the flocculent precipitate of ‘humic 
acid’ forms. The mixture is filtered into a 50 ml. volumetric flask, and 
the precipitate is washed with 0.01 N hydrochloric acid until the volume of 
filtrate is about 40 ml. The filtrate is adjusted to pH 3.0, and diluted to 
volume after the addition of 2 ml. of acid molybdate reagent. 


D. Colorimetric Procedure 


The photoelectric colorimetry of phosphorus in soil extracts has been dis- 
cussed at length elsewhere (4, 9), and only an outline of the necessary manipu- 
lations is given here. 

A 10 ml. aliquot of the final solution, containing acid molybdate reagent, 
is pipetted into a colorimeter tube. The tube is placed in the colorimeter, in 
which the 660 my filter has been inserted, and the rheostat so adjusted that 
the galvanometer shows a reading of 100*. One drop of 1% stannous chloride 
reagent is then added and mixed with the contents of the tube, and the gal- 
vanometer reading is recorded when the colour has developed to a maximum, 
which is usually about 10 min. after the addition of stannous chloride. The 
minimum galvanometer reading, corresponding to maximum colour, gives 
directly the concentration of phosphate phosphorus at the final dilution, by 
reference to calibration data. The appropriate calibration constant is 
obtained by measuring the light absorption in a phosphate solution of known 
concentration which contains the same concentration of ammonium chloride 
as the unknown. 

The results in this paper are expressed as milligrams of phosphorus per 
gram of soil. 


* In reading a series of tubes it is convenient at this point to record the ‘‘centre setting”, 1.e., 
the galvanometer reading when the tube is removed. After the colours have been developed, adjust- 
ment to the appropriate centre setting is made for each tube before it is placed in the colorimeter 
for the final reading. 
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Results and Discussion 


Some results obtained by the use of the method are presented in Table I. 
It may be seen from these data that a considerable proportion of the organic 
phosphorus is brought into solution by the acid extraction. On the other 
hand, following the action of relatively strong acid, the ammonia solution dis- 
solves only slight traces of inorganic phosphorus. The sum of the amounts 
of phosphorus dissolved by the two extractions approaches the total phos- 
phorus of the soil, the percentage extraction being from 90 to 100% in the 
typical instances reported. The almost complete recovery of the phosphorus 
in the muck soils indicates that the extraction of organic phosphorus is com- 
plete. The agreement between duplicate values for the proportion of the 
soil phosphorus that is organically combined appears to be very satisfactory, 
when it is considered that these values may combine the errors of four colori- 
metric phosphorus determinations. The analysis of an ignited sample offers 
convincing proof that the method gives a true representation of the organic- 
inorganic phosphorus distribution. 

Further extraction with ammonia of the residue from mineral soils brings 
no more phosphorus into solution. It may be assumed that the phosphorus 
not dissolved by the two extractions is inorganic phosphorus contained in 
undecomposed mineral fragments. In one case the phosphorus content of 
the grey, silicious residue from the extraction of a mineral soil was determined, 
after magnesium nitrate fusion, and proved to be about 6% of the total soil 
phosphorus. This result isin harmony with the data for percentage extraction. 

In order to ensure that no interfering substances or conditions were affecting 
the direct determinations of inorganic phosphorus, known amounts of phos- 
phate were added to some of the extracts, and the recovery of these known 
additions was noted. It was found that the recovery was 99 to 100% in the 
acid extracts, but rather less satisfactory in the alkaline extracts, ranging 
from 95 to 100%. It was thought that the “humic acid” precipitate might 
be adsorbing inorganic phosphate. Some of these precipitates were dissolved 
in dilute ammonium hydroxide, reprecipitated, and the phosphate in the 
filtrate was estimated. The merest traces of phosphate were found, and it 
was concluded that inorganic phosphate is not adsorbed to an appreciable 
extent. Inorganic phosphate was determined directly in a dilute solution 
of the alkaline extract in some instances where the dissolved organic matter 
did not give too flocculent a precipitate upon acidification of the solution. 
The results were almost identical with those given by the filtration method. 
The difficulty in obtaining complete recovery of the phosphate in the alkaline 
extracts appears to be due mainly to the high concentrations of ammonium 
chloride present in these solutions. As shown elsewhere, ammonium chloride 
in high concentration not only depresses the colour intensity but also renders 
the colour less stable, so that it is almost impossible to perform quantitative 
colorimetry. However, the amount of phosphate involved in this particular 
determination is such a small proportion of the total phosphorus that the 
effect upon the result of the method as a whole is negligible. 
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It is known that organic phosphorus compounds are fairly stable in both 
acid and alkaline solutions at ordinary temperatures (7), and Schollenberger (8) 
showed that dilute ammonium hydroxide does not have any destructive action 
on the organic phosphorus compounds in soil. It was nevertheless thought 
important to test the stability of soil organic phosphorus in the presence of the 
relatively strong acid used for extraction. A series of analyses of an acid 
extract showed no change in the phosphorus distribution in a period of 10 days. 
The acid concentration of the final extract is, however, only 1 N, and an experi- 
ment was designed to compare the effects of different concentrations of acid 
applied directly to the soil in the usual procedure. 

Hydrochloric acid was used at the following concentrations: 0.5 N, 1 N, 
3.N,5 N. The results are summarized in Table II. 

There does not appear to be any consistent difference between the results 
obtained with different concentrations of acid. The use of 4 N acid seems to 
be quite reasonable in the light of these data; in fact, it is evident that the 
concentrations of the extracting acid may vary between wide limits without 
affecting the determination. It is felt that this fact is in strong support of the 
validity of the present method. 

Table II shows that the total phosphorus extracted is appreciably more 
when 3 N or 5 N acid is used. The solubility of organic phosphorus in the 
acid solution increases progressively with the concentration of acid, and 
ammonia extraction removes the remainder in each case. The amount of 


TABLE II 


INFLUENCE OF THE CONCENTRATION OF THE EXTRACTING ACID ON THE DETERMINATION 
OF TOTAL ORGANIC PHOSPHORUS 


li 
Ceinenieaiinn Milligrams of phosphorus per gram of soil Org. P, 
of extracting Per cent} % of 
acid Acid extract Ammonia extract Total Total pong total 
(HC) P org. P 
Total | Inorg. Org. Total | Inorg. Org. ext’d Pp 
O.5N 0.255 0.207 0.048 0.455 0.052 0.403 0.710 0.451 91.0 57.8 
1 N 0.294 0.234 0.060 0.415 0.044 0.369 0.709 0.429 90.9 55.0 
3 0.467 | 0.287} 0.180 | 0.290 | 0.020| 0.270| 0.757 | 0.450| 97.0 57.7 
5 0.555 | 0.322} 0.233 | 0.205 | 0.012 | 0.193 | 0.760] 0.426] 97.4 54.7 


* Soil No. 8, a podsol, with a total P content (air dry) of 0.078%. 


inorganic phosphate in the ammonia extract decreases as the concentration 
of the extracting acid is increased. This observation favours the extraction 
with strong acid (3 to 5 N) since its use minimizes the errors of colorimetry 
in regard to the alkaline extract. Indeed, when the acid used is 3 N or 
stronger, only a minor error would be incurred by assuming that all the phos- 
phorus of the ammonia extract is organic phosphorus. This suggests a means 
by which a very good estimate of the organic phosphorus may be obtained 
without the use of a photoelectric colorimeter, as visual colorimetry is feasible 
in the acid extract. 
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The results in Table I suggest an interesting relation between soil type and 
the proportion of phosphorus organically combined. In muck soils, which 
owe their existence to the inhibition of organic decay, the proportion of organic 
phosphorus is, very high. In the podsols, where the organic matter is poorly 
decomposed, the proportion of organic phosphorus is intermediate, and in the 
brown forest soils, where the decomposition of organic matter is well advanced, 
the proportion of organic phosphorus is lowest. It is, perhaps, of interest 
that soil No. 6 is a very fertile soil, probably the most fertile of the group 
analyzed. However, the data at hand are too few to permit any definite 
conclusions to be drawn in this connection. 
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PRESSURE, VOLUME, TEMPERATURE RELATIONS OF 
ETHYLENE IN THE CRITICAL REGION. I.! 


By J. Dacey’, R. McINTosH’, AND O. Maass* 


Abstract 


Pressure isothermals of the system were determined both above and below the 
temperature of disappearance of the meniscus. Regions of constant pressure 
with changing velume were observed with all isothermals, in agreement with the 
prediction of Mayer and Harrison (2, 3). The complete envelope of this region 
above the critical temperature has not as yet been determined. The pressure 
of the heterogeneous system at 9.50° C. was found to correspond within experi- 
mental error to the pressure of the flat portion of the homogeneous isothermal at 
the same temperature. The pressure of this heterogeneous system was found to 
be independent of the mass-volume ratio; this is in agreement with the experi- 
mental behaviour previously discovered by the authors (5) for such a system at 
9.80°C. Moreover, the pressure of the heterogeneous system was found to be 
identical with that of the system which had been heated at constant volume to 
a temperature at which the density difference in the tube is known to have been 
destroyed, and then cooled to the same temperature of 9.50° C. 

The isothermal at 9.60°C. was studied, and shown to be the equilibrium 
curve at that temperature. A hysteresis was observed on reversing the direction 
of measurement, 1.e., from “‘vapour”’ to compressed “‘liquid”’. This hysteresis was 
found to be caused by the time lags which are observed in passing from the 

“vapour” region of an isothermal to regions of high density. Certain apparent 
discrepancies between the work of the present authors and that of Geddes and 
Maass have been observed and cannot as yet be explained. 

The phenomena observed are interpreted on the basis of a difference sisiaaiin 
the gaseous and liquid states of aggregation, with a structure assigned to the 
latter. 


Introduction 


The present work is a continuation of that reported by McIntosh and 
Maass (5) on the density, pressure, and temperature relations of ethylene, 
and has been directed primarily toward the study of certain predictions 
made by Mayer and Harrison (2, 3) of the behaviour of condensing systems, 
and of the hypothesis of Maass that the liquid state of aggregation may persist 
above the temperature of disappearance of the meniscus and is characterized 
by a dynamic regional orientation. ~- 


Specifically, Mayer and Harrison have predicted the existence of a region 


above the critical temperature throughout which (3), = 0. This region 


has an upper limit at that isothermal for which (37), = 0 for one point 


only. The isothermals of Geddes and Maass (4) showed no such region. 
This can be explained on the assumption that the experimental points of the 
isobars, from which the isothermals were evaluated, were not taken with 
sufficiently small temperature increments. 


1 Manuscript received February 28, 1939. 
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In the present work, isothermals were determined experimentally at volume 
intervals as minute as possible, particularly throughout the volume ranges 
where flat portions of the curves might be expected. Relative pressures are 
determined to the limit of error that is inherent in the manual thermostatic 
control of the temperature of the carbon dioxide, that is, 0.01 atm. Within 
this limit the prediction of Mayer and Harrison has been verified. The data 
of these and other experiments are given below. 


Experimental 


The apparatus used in this investigation has been adapted from the original 
design of Geddes and Maass (4). It consists of two Pyrex glass bombs, 
supported on brass rods, and connected from the bottoms by a U-tube of 
heavy walled Pyrex capillary. The apparatus is filled with pure mercury 
' to about the half-way position in each bomb. This mercury acts as a piston 
in the manipulation of the ethylene, which is confined in the right-hand bomb, 
and also transmits the pressure generated by the carbon dioxide in the left- 
hand bomb. The carbon dioxide is added in the proper amount so that a 
two-phase system will be maintained through the volume and temperature 
ranges required. In these respects the apparatus of Geddes and Maass has 
been copied exactly. In the present arrangement, however, the spiral and 
float have been omitted, and differences in volume are recorded by the position 
of the mercury level as read upon a cathetometer scale. The internal diameters 
of the tubes were calibrated by measuring the height of known weights of pure 
mercury, and were found to be constant within the experimental error of the 
calibration. The determination of the absolute volume necessitated a 
geometric estimation of the volume of the draw-off; this was checked by 
repeating the determination of a point upon the isobars of McIntosh and 
Maass (5). The agreement was very good. 

The absence of the spiral and float permitted the use of much narrower 
tubing, with consequent decrease in the danger of explosions. The inability 
with this design, to detect the presence of heterogeneous systems is not a 
great disadvantage, because the conditions under which these systems exist 
were well known from previous work, and the pressures of such systems were 
the point of interest rather than the actual density differences. Moreover, 
it was now possible to stir the ethylene by means of a magnetic stirring 
arrangement. 


Experimental Results 


Five isothermals of the system ethylene were obtained and are represented 
in Fig. 1 and the values tabulated in Table I. The experimental procedure, 
the same in all cases, was as follows. The ethylene was compressed to a high 
density at 8.9° C. by adjustment of the carbon dioxide pressure. The tem- 
perature of the ethylene was then raised at constant volume to the temper- 
ature of the desired isothermal by balancing the change in pressure by further 
careful adjustment of the carbon dioxide pressure. After equilibrium had 
been attained, small decreases in the pressure were made and the equilibrium 


PRESSURE-VOLUME-TEMPERATURE DATA 


TABLE 
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Volume, Pressure, 
ee. atm. 


Volume, 
cc. 


Pressure, 
atm. 


Volume, Pressure, 
CC. atm. 


Temperature, 8 .92° C. 


Temperature, 9.22° C. 


Temperature, 9.42° C. 


2.532 49 .93 
2.587 49.73 
2.672 49.54 
2.893 49.45 
3.271 49.45 
3.681 49.46 
4.303 49 
4.317 49.34 


NNN 
oo 
w 


2.710 50.26 
2.735 50.15 
2.838 50.06 
2.950 50.02 
3.176 49.99 
3.300 49 .99 
3.550 49.97 
3.123 49.94 
3.910 49.90 
4.104 49.85 


Temperature, 9.50° C. 


Temperature, 9.60° C. 


2.705 50.42 
2.122 50.32 
2.762 50.22 
2.789 50.19 
2.829 50.17 
2.908 50.12 
2.970 50.10 
3.235 50.08 
3.455 50.07 
3.624 50.06 
3.857 50.03 
4.197 49.90 


w 


4.000 50.10* 
3.958 50.16 
3.903 50.20 
3.782 50.23 
3.586 50.29 
3.148 50.34 
2.858 50.36 
2.753 50.42 


* Increasing pressure, thermal equilibrium only established. 
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49.3 


8.92°C. 


25 


3.0 3.5 
VOLUME, CC. 


Fic. 1. 


: 
208 
50.18 
50.06 
49.96 
49.87 
49.85 
49.80 
49.75 
49.75 
49.73 
49.72 
4.006 49.68 
4.174 49.65 
Temperature, 9.60° C. | 
13 50.45 
56 50.34 
54 50.27 
35 50.18 
it 60 50.17 
pat 12 50.15 
07 50.14 
10 50.13 
59 50.09 
50.5) 
960°C. 
9.50°C. 
| | | 942°C. 
922°C. 
| | 40 
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volumes determined from the mercury level as read upon the cathetometer 
scale. In that region where a horizontal part of the curve is shown, the 
decrease in pressure is exactly compensated for by the change in mercury 
level. 


All isothermals showed a region wherein (37) = 0, although the 
T 


isothermals at 9.50°C. and 9.60°C. are respectively at and 0.10° higher 
than the temperature at which the meniscus disappears. The envelope of 
this region has not yet been completely determined, but insofar as the work 
has been carried out, the prediction of Mayer and Harrison mentioned above 
has been verified. 


Isothermal at 9.60° C. 
Since the deduction of the existence of a two-phase system from the observa- 


tion of a region wherein (37) = Ois open to the objection that the behaviour 
T 


of a one-phase system, if the compressibility is great enough, may simulate 
that of a two-phase system, and since the compressibility of liquids near the 
critical point is admittedly enormous, the direction of measurement of the 
points along the isothermal at 9.60° C. was reversed. That is, when the 
substance had been expanded to the last point on the isothermal, the pressure 
was increased. As soon as temperature equilibrium had been established, the 
pressure was maintained constant for 10 min. By following this procedure, 
the upper curve of Fig. 2 was obtained. This showed the compressibility 


PRESSURE , ATM, 


35 
VOLUME , CC. 
Fic. 2 


of the system to be dependent upon its previous history. That the points 
of the lower curve are not dependent on the previous condition of the system, 
provided only that expansion into the true “‘gaseous” region has not occurred, 
was shown by repeating, several days later, the determination of a point on 
the flat portion of the curve. The procedure was as follows. The system was 
compressed to a high density, the temperature raised as previously described, 
and the pressure adjusted to the value corresponding to the flat portion of 
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the curve. Immediate expansion took place, and an equilibrium value was 
attained at the volume represented by the shaded circle. This point showed 
no change after more than 15 min. under these conditions. The pressure 
was then adjusted in order to repeat the last ‘‘vapour’’ point shown upon the 
diagram. This point (also shaded circle) was exactly repeated, and no 
tendency to change was exhibited at the end of 15 min. 

The question whether the pressure on the higher curve would fall spon- 
taneously to the value on the lower curve was then investigated. After 
expansion to the “‘gaseous”’ region of density the pressure on the system was 
adjusted to give the triangular point on the upper curve. When this point 
had been reached, the volume was maintained as constant as possible, and a 
pressure-time curve obtained. It was found that the pressure fell quite 
rapidly to the value of the original curve, where it remained constant. The 
difficulty of keeping the volume constant is very great, as the limit of adjust- 
ment of pressure is 0-01 atm., which represents a change of nearly 4 mm. in 
the height of the mercury. However, the volume can be kept constant 
enough to show clearly the true trend of the pressure with time. The figures 
are given in Table II. 


TABLE II 
Volume Temperature Pressure Time 
3.59 9.60 50.28 0 
3.55 9.60 50.27 7 
3.55 9.60 50.24 13 
3.55 9.60 50.22 20 
3.57 9.60 50.20 26 
3.5 9.60 50.18 34 
9.60 50.18 40 


Note: Zero time was taken at about 12 min. after the 
original pressure adjustment was made. 


The information obtained from these experiments indicates the interesting 
deduction that, as shown by Geddes and Maass, although compression 
destroys the macroscopic heterogeneity of the system, the two-phase character 
is re-formed by compression, even at a temperature 0.10° higher than the 
temperature of disappearance of the meniscus. This leads to the conclusion 
that over short ranges of temperature above that at which the meniscus 
becomes invisible, the nature of the system is determined by the density, 
or the mass-volume ratio. Further support of this statement may be derived 
from the following experiments. 


Pressure Measurements of the Heterogeneous and Homogeneous Systems at 9.50° C. 

The question of relative stability of the heating and cooling curves of such 
a typical isochore as appears in Fig. 3, as measured by the pressure of the 
systems over a range of density from below the critical density to well above 
it, was determined by the following series of experiments. The volume of 
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the system was fixed by freezing the mercury in the centre U-tube at a value 
of the average density where a heterogeneous system is known to exist from 
the previous work of McIntosh and Maass. The temperature of the ethylene 
was then raised to 9.50° C., the mercury melted, and the pressure of the system 
corresponding to a point upon the heating curve of the isochore determined. 


= 9 10 i 12 


TEMPERATURE , °C, 


Fic. 3. 


The mercury was then frozen, and the temperature raised to 11.5° C., where 
the density of the system has been shown to be constant. After sufficient 
time had been allowed for the medium to reach equilibrium, the temperature 
was again lowered to 9.50° C. and the pressure determined. It was shown 
that over the range of densities studied, the pressures shown by the upper 
curve, A, and lower curve, B, are identical, and are the same as the pressure 
corresponding to the 9.50° C. isothermal of the homogeneous system. The 
data are given in Table III. The pressures, although recorded to only 0.01 


TABLE III 
Av. Dens. Pressure, Pressure, Temp. 
. 1990 50.08 50.08 9.50 
° . 2040 50.08 50.08 9.50 
. 2090 50.08 50.08 9.50 
.2173 50.08 50.08 9.50 
. 2287 50.08 50.08 9.50 


Pressure of flat portion of homogeneous isothermal at 
9.50 is 50.07 atm. Range of density of horizontal portion 
is 0.229 to 0.198 gm. per cc., approximately. 


atm., are more nearly equal than 0.01 atm. The reason for this statement 
is that, although changes in pressure can be determined with any degree of 
certainty only to that limit, smaller differences in pressure are shown by 
the change in mercury level on the cathetometer scale. That is, when a 
pressure measurement had been made during the determination of the heating 
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curve, A, the mercury frozen, and the conditions adjusted to give the point 
upon the cooling curve, B, and the mercury melted, no change in mercury 
level could be detected by means of the cathetometer. 


The agreement between the pressure of the macroscopically heterogeneous 
system and that corresponding to the horizontal portion of the isothermal 
of the homogeneous system, together with the further similarity of density 
range throughout which the heterogeneous system is known to exist on the 
one hand, and throughout which the pressure of the homogeneous system is 
known to remain constant on the other hand, indicates that for these tem- 
peratures which exhibit horizontal regions on the pressure-volume curves, 
heterogeneous systems may be obtained, which are apparently stable and 
will show no tendency to change spontaneously to systems of a homogeneous 
character. 


Further, apparently then, the flat portion of the cooling curve of a typical 
isochore cannot be regarded as representing a homogeneous gas system with 
the liquid structure destroyed, for such a gas system would have a pressure 
independent of the volume over quite appreciable ranges. It is concluded, 
therefore, that the system, although macroscopically homogeneous, is still a 
two-phase system, with the effect of the temperature being to disperse, or 
mix mechanically, the two phases. This opinion is in agreement with the 
conclusion reached by Naldrett, Mason and Maass from a study of critical 
opalesence which is to be published shortly. It also explains the observation 
of McIntosh and Maass (5) that the phase densities of ethylene at the critical 
temperature are dependent on the mass-volume ratio of the system. From 
this new point of view, the phase densities will be dependent on the amount 
of dilution of one phase by the other, which, as the meniscus rises or falls 
rapidly near the critical temperature, depending upon the value of the average 
density, will be a function of the mass-volume ratio. 


Effect of Mechanical Stirring upon the Pressure of the Heterogeneous. System. 


In the design of apparatus used in this investigation it was possible to include. 
a magnetic stirring arrangement in the ethylene side of the system. This was 
done in order to ascertain whether mechanical stirring, which has been 
emphasized as necessary in criticisms of work done in this laboratory, would 
affect the phenomena which are observed. 


A heterogeneous system was therefore obtained at 9.50° C. and its pressure 
recorded. This system was then violently stirred by movement of the 
magnet up and down the length of the tube for a period of 20 sec. Several 
minutes were allowed to elapse, and the pressure was again recorded. This 
treatment was repeated four times. No change in pressure could be detected 
and the mercury level as read upon the cathetometer scale was unchanged. 
Although no mass movement of material could be detected during the stirring, 
the stirring was of sufficient violence to cause large visible transfers of liquid 
short distances into the vapour phase when the temperature was below the 
critical. This experiment lends strong support to the statement of McIntosh 
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and Maass that “mechanical stirring that does not cause local heating and 
cooling, or expansion and compression, does not influence the phenomena 
observed, either in kind or in rate of attainment of equilibrium.” 

Moreover, ‘the fact that the phenomena reported in this paper, which are 
similar to those reported in previous publications from this laboratory, was 
observed in a tube that contained no spiral, completely refutes the argument 
of Ruedy (4) that the presence of a spiral acts as a deterrent to the proper 
mixing and diffusing of the system, and is the reason for the anomalous 
observations that have been made in the critical region. 


Apparent Discrepancies between Work of Geddes and Maass and Present Work 


A disagreement in the values of the critical pressure of ethylene has been 
found. Geddes and Maass found that the critical pressure depended upon 
the mass-volume ratio and reported the two values of 49.98 and 50.01 atm. 
at the average densities of 0.2056 and 0.2302 gm. per cc. respectively (5). 
The present work gives a critical pressure of 50.07 atm. over a density range 
of 0.198 to 0.220 gm. per cc. approximately. No attempt to explain this 
difference is considered feasible at the present stage of the work. 
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